Abstract. Long-slit echellograms of the high excitation 
planetary nebula NGC 1501, reduced according to the 
methodology developed by Sabbadin et al. (2000a, b), al- 
lowed us to obtain the "true" distribution of the ionized 
gas in the eight nebular slices covered by the spectroscopic 
slit. A 3-D rendering procedure is described and applied, 
which assembles the tomographic maps and rebuilds the 
spatial structure. The images of NGC 1501, as seen in 12 
directions separated by 15°, form a series of stereoscopic 
- - - pairs giving surprising 3-D views in as many directions. 
. The main nebula consists of an almost oblate ellipsoid of 
■ moderate ellipticity (a~44 arcsec, a/b~1.02, a/c~l.ll), 
fS| brighter in the equatorial belt, deformed by several bumps, 
and embedded in a quite homogeneous, inwards extended 
q_) cocoon. Some reliability tests are applied to the rebuilt 
pLn ■ nebula; the radial matter profile, the small scale density 
. fluctuations and the 2-D (morphology) - 3-D (structure) 
p\J ' correlation are presented and analysed. The wide applica- 
, tions of the 3-D reconstruction to the morphology, phys- 
' ical conditions, ionization parameters and evolutionary 
status of expanding nebulae in general (planetary nebulae, 
nova and supernova remnants, shells around Population I 
^f) Wolf-Rayet stars, nebulae ejected by symbiotic stars, bub- 
C^l bles surrounding early spectral type main sequence stars 
,; etc.) are introduced. 
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1. Introduction 

The planetary nebula (PN) phase (common to most of 
intermediate and low mass stars in their late evolution; 
Iben, 1984; Vassiliadis & Wood, 1993; Schonberner, 1997) 
is characterized by an extreme variety of morphologies: 
disks, butterflies, point-symmetric or multiple shells, hour- 



glasses etc. ( 


Greig, 1972|; 


Bahck, 1987 




Schwarz ct al., 


1992; 


Manchado et al., 1996). 







Different mechanisms, acting before, during and after 
the nebula ejection, can contribute to generating the ob- 
served shapes. These include the mass-loss history, stel- 
lar rotation and/or precession, magnetic fields, duplic- 
ity of the central star, wind interactions, hydrodynamical 
processes, ionization and recombination (Tylenda, 1986; 
Morris, 1987|; [Pascoh, 1992| ; [Pwarkadas fc Balick, 19981 ; 
Garcia-Segura et al., 1999| ). 

In order to disentangle the physical processes forming 
and shaping a PN, the first - and most difficult - observa- 
tional step is the nebula de-projection, i. e. the recovery of 
the ionized gas 3-D structure. In this paper, which is the 
result of a long-standing effort (Sabbadin, 1984, Sabbadin 
et al, 1985, 1987, 2000a, 2000b), we obtain the spatial 
distribution of the PN matter by means of high resolution 
spectroscopy. 

As a test case we selected NGC 1501 for which we have 
already presented a partial tomographic analysis (Sab- 
badin et al. 2000b, hereafter PI). In PI, long-slit echel- 
lograms at four position angles gave the radial density pro- 
file, the mean electron temperature and turbulence and a 
rough 3-D model. The accurate spatial reconstruction was 
deferred to a more complete spectroscopic coverage. 

Recently we have observed NGC 1501 at four addi- 
tional position angles and developed a rendering proce- 
dure which integrates the tomographic slices and rebuilds 
the detailed spatial distribution of the ionized gas. 

This work, in many respect the logical continuation 
of PI, is structured as follows: the nebula and the ob- 
servational material are introduced in Sect. 2; in Sect. 3 
we illustrate the reduction procedure, in particular, how 
to derive the electron density profile from the Ha flux of 
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Fig. 1. Spectroscopic slit positions and position angles, su- 
perimposed on a broad-band R frame of NGC 1501 taken 
with the 3.58m Itahan National Telescope (TNG). The 
seeing was FWHM=0.76 arcsec. North is up and East is 
to the left. 



the zero- velocity pixel column; Sect. 4 contains the to- 
mographic maps and Sect. 5 describes and analyses the 
spatial model; a general discussion is presented in Sect. 6 
and conclusions are drawn in Sect. 7. 



2. The nebula and the observational material 

The high excitation PN NGC 1501 (PNG144.5-)-06.5, 
Acker et al., 1992) is "a mass of irregu lar condensations, 
with some big holes" ( Minkowski, 1968 ) " bearing a res em- 
blance to the convolutions of the brain" ( Pease, 1917 ). 

The literature-based data for the whole system (nebula 
-|- exciting star) are summarized in Table 1. 

AA4500-8000 A, spatially resolved spectra of NGC 1501 
(including flat field, Th-Ar calibration and spectropho- 
tometric standard) were secured with the Echelle spec- 
trograph attached to the Cassegrain focus of the 1.82m 
Asiago (Italy) telescope and equipped with a Thompson 
1024x1024 pixels CCD (spatial scale=0.79 arcsec pix~^; 
spectral scale=9.0 km s~^ pix"-'^). 

Eight position angles were covered in two observing 
runs (the slit positions are shown in Fig. 1 superimposed 
on a broad-band R frame of the nebula): 

Run 1) P.A. = 10°, 55°, 100° and 145° on Decem- 
ber 1998 (non photometric nights; seeing 1.4 - 1.7 arcsec; 
spectral resolution 13.5 km s^^; spectra used in PI); 

Run 2) P.A. = 30°, 80°, 125° and 170° on December 
1999 (photometric nights; seeing 1.2 - 1.4 arcsec; spectral 
resolution 12.0 km s""'^). 

In all cases the exposure time was 1800 seconds. 

A detailed description of the spectroscopic character- 
istics of NGC 1501 is given in PI. In short: it is a high 
ionization, density bounded PN; stratification effects are 
negligible; the position-velocity maps can be fitted by a 
triaxial ellipsoid deformed by several hemispheric bubbles 
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coordinates (2000.0) 
apparent size [arcsec] 
log F(H^) [mWxm-2] 
c(H/3) 

excitation class 
individual distance [Kpc] 

statistical distance [Kpc] 



adopted distance [Kpc] 
mean linear radius [pc] 
expansion velocity [km s~^] 

turbulence [km s~^] 
structure 



NEBULA 
a: Q4h 06"" 59.7" 5: +60" 55' 14" 
56x48 (Curtis, 1918); 68 faint, roundish halo (PI) 
-11.28 (Collins et al., 1961); -11.20 (Acker et al., 1991) 

0.96 (Kalcr, 1976); 1.1 (Stanghellini et al., 1994); 1.11 (Ciardullo et al., 1999); 1.05 (PI) 
8 pecuhar (Page, 1942) 

2.0 (Acker, 1978); 1.4 (Pottasch, 1983); 1.2 (Sabbadin, 1986) 
(all based on the average extinction in the galactic disk) 

1.45 (O'Dell, 1962); 1.20 (Cahn & Kaler, 1971); 1.78 (Cudworth, 1974); 1.16 (Cahn, 1976); 

1.16 (Acker, 1978); 1.10 (Daub, 1982); 1.10 (Maciel, 1984); 0.90 (Ammicl ct al., 1984); 

1.50 (Sabbadin, 1986); 1.16 (Cahn et al., 1992); 1.31 (Van de Steene & Zijlstra, 1994); 

1.21-1.75 (Zhang, 1995) 

1.30 

0.164 

39 in [QUI] (Robinson et al., 1982); 38 in [OIII] and HI (Sabbadin & Hamzaoglu, 1982); 

40 in Hell, [OIII], HI and [Nil] (Neiner et al., 2000); 38 to 55 in [OIII] and HI (PI) 
10 (Neiner et al., 2000); 18 (PI) 

prolate spheroid of moderate ellipticity (Sabbadin & Hamzaoglu, 1982); 

thin ellipsoid deformed by a pair of large lobes along both the major and the intermediate 

axes and by a heap of bumps spread over the whole nebular surface (PI) 



mv 

spectrum 

logTstar 

mass loss rate [Moyr~^] 
terminal wind velocities [km s~ 
note 



CENTRAL STAR 

14.23 (Tylondact al., 1991); 14.36 (Ciardullo ct al., 1999) 

WC6 (Swings & Swensson, 1942); WC-OVI (AUer, 1976); WC4 (Tylenda et al., 1993) 
4.98 (TZife//; Sabbadin, 1986); 4.91 (TZne/j; Stanghellini et al, 1994); 
5.13 (atmosphere model; Koesterke & Hamann, 1997) 
5.2x10"'' (Koesterke & Hamann, 1997) 

1800 (Koesterke & Hamann, 1997); 3260 to 3460 (Fcibclman, 1998) 

the star is a nonradial g-mode pulsator with periods ranging from 5200s to 1154s 

(Bond et al., 1998; Ciardullo & Bond, 1998) 



Further informations on NGC 1501 are contained in Pottasch (1983) and Acker et al. (1992). 



protruding from the central figure; the main emission oc- 
curs in a narrow shell whose density profile is steep out- 
wards and flatter inwards. 

We believe that a general observational note is appro- 
priate here, concerning the common employment by most 
authors of an interference filter when obtaining echelle 
spectra of PNe (to isolate a single order containing one or 
more interesting emissions, like Ha + [Nil], and to avoid 
the superposition of orders). No filter was used in our ob- 
servations, for the following reasons: 

- PNe present a discrete spectrum and the superposition 

of lines belonging to different orders is quite improba- 
ble (even with a slit covering several orders); 

- the introduction of an interference filter, reducing the 

spectral range from several hundred A to a few dozen 
A, drastically limits the number of usable emissions 
per frame; 

- a careful reduction procedure allows us to take into ac- 

count all the possible sources of inaccuracy (line cur- 
vature, order tilting, optical distorsions, spectral sen- 
sitivity, CCD cosmetics) over the whole frame, thus 
obtaining the flux and wavelength calibrated, spatially 



resolved structure of all emissions present in the echel- 

logram (some dozens in the case of the brightest PNe, 
representing the main ionic species). 



3. Reduction procedure: the radial electron 

density (ATe) profile from the Ha absolute flux 
in the zero-velocity pixel column 

The spectral reduction process was the same as described 
in PI, but for the absolute calibration of the radial iVe 
proflle. Due to night sky variability, in PI we were forced 
to obtain the electron density of NGC 1501 by comparison 
with NGC 40 (both nebulae observed on the same nights 
and with the same instrumental setup). This was unnec- 
essary for "run 2" , since the sky was photometric and the 
proper nebular flux calibration was performed using the 
spectrophotometric standard. 

Thus, the starting point is the classical expression 

rR 

ATTD'^F{Ha) = hv3,2 a3,2 N{H+)Nee-iTrr^dr (1) 

^0 

where: 
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Fig. 2. Radial electron density profiles at the four position 
angles of "run 2" , obtained from the absolute flux in the 
Ha zero-velocity pixel columns. 

D is the distance, F{Ha) is the Ha flux corrected for 
extinction, is the Ha frequency, as. 2 is the Ha re- 
combination coefficient, N{H^) is the ionized hydrogen 
density and e is the filling factor, representing the fraction 
of the nebular volume which is filled by matter at density 
Ne. 

We restrict the integration to the portion of the nebula 
identified by the zero- velocity pixel column (corresponding 
to the gas which is expanding perpendicularly to the line 
of sight). 

The Ha flux (Fj) in each pixel of the zero- velocity col- 
umn, corrected for interstellar absorption and instrument 
resolution (for details see Sabbadin et al. 2000a), is given 
by: 



A-kD^Fi = ATTiHo,N{H+)NeViti 
where: 



[ergs 



(2) 



jnce is the Ha emission coefficient for the case B of 
Baker fc Menzel (1938| ); for typical nebular densities 
its dependence on Te is obtained by interpolation of 
Brocklchurst (1971 ) data: AirjHa = hv^ 2 = 0.20 x 

10-20yg-0.94. 

Vi is the "local nebular volume" defined as V; = Ax Ar, 
A being the pixel area and Ar — R x (AV/Vcxp) the 
radial thickness of the zero-velocity pixel column; 

ei is the "local filling factor" , indicating the fraction of 
the local nebular volume, Vi, which is actually filled by 
matter with density A'^e. 



The radial electron density profiles of NGC 1501 de- 
rived from eq. (2) at the four position angles of "run 2" 
are shown in Fig. 2 (assuming Te = 11500°K, Ne ~ 1.15 
iV(H+), c(H/3) = 1.05, D = 1.30 Kpc and, in first approx- 
imation, e;=l). Their complex structure confirms (at least 
qualitatively) the results already obtained in PI; in par- 
ticular, the main nebula consists of a thin shell presenting 
large inwards tails whose extent is anti-correlated to the 
height of the density peak. 

A quantitative Neruni—Nerun2 comparison can be per- 
formed only at the intersection of the eight slit positions, 
i. e. at the apparent stellar position; this needs the ac- 
curate tomographic reconstruction illustrated in the next 
Section. 

4. Tomography 

The tomographic analysis of PNe, originally introduced 



and applied to plate spectra by Sabbadin (1984) and Sab 



badin et al. (1985, 1987), has been recently extended and 
refined by Sabbadin et al. (2000a, b). Its application to 
the Ha emission of "run 2" spectra allows us to derive the 
electron density maps shown in Figure 3. Following PI, 



we adopted a linear distance- velocity relation, an electron 
temperature of 11500°K, a turbulence of 18 km s~"'^ and a 
radius of the shell peak at the position of the central star 
(in radial direction) of 20 arcsec. 

The maps of Fig. 3, combined with the similar maps of 
Fig. 5 of PI, give the spatial iVe distribution in the eight 
nebular slices covered by the spectroscopic slit. NGC 1501 
being a high ionization, density bounded PN, the A^(H+) 
and A^(0++) tomographic maps coincide with the Ne ones 
(but for a scaling factor of 1.15 and 3.8 x 10^, respectively; 
see PI). 

The comparison of the electron densities in the neb- 
ular regions common to the eight position angles (i. e. 
at the apparent stellar position) gives a very satisfactory 
agreement (to within ±3%) between Nemni and A^erun2, 
confirming the validity of the "escamotage procedure" 
adopted in PI. Our satisfaction is mitigated by the consid- 
eration that, in any case, the inaccuracy in the absolute 
A^e calibration based on eq. (2) is of the order of 10%, 
mainly due to the uncertainty in the distance of the neb- 
ula. 

The only direct Ne determination in NGC 1501 us- 
ing diagnostic line ratios dates from AUer & Epps (1976). 
These authors analysed at low spectral resolution a small 
region located at a distance of 20 arcsec from the central 
star (in PA= 125°) and derived /(6717)//(6731) = 0.78, 
corresponding to Ne = 1200 cm'^ (for Te = 11500°K). 
From the tomographic map at PA= 125° (Fig. 3) we ob- 
tain, at a distance of 20 arcsec in the E-SE sector, density 
peaks of 950-970 cm'^, i.e. our Ne{SB) {SB stands for 
"surface brightness" ) is 20% lower than the corresponding 
[SIX] electron density reported by AUer & Epps (1976). 



Possible causes of this discrepancy are (in order of im- 
portance): 

a: small-scale density fluctuations. We adopted a local 

filling factor ei — 1; actually, the surface brightness 

1/2 

method (cf. eq. (2)) gives A^e x . On the other 
hand, the [SIX] line intensity ratio measures the true 
electron density of the brightest (i. e. densest) nebular 
regions; thus the iVe(5'i?)/A^e([SXX]) ratio would indi- 
cate a local filling factor ei ~ 0.66 in the main shell of 
NGC 1501; 

b: distance inaccuracy. From eq. (2) we have that 
Ne(SB) oc {l/ny/^; to obtain Ne{SB) = A^e([SXX]) 
we need £>ngci501 — 0.83 Kpc (assuming e; — 1). This 
value is lower than any of the estimates reported in 
the literature (see Table 1); 

c: ionization. NGC 1501 being a high ionization PN, S+ 
is localized in the outer parts of the densest regions, 
where the shadowing of the UV stellar flux decreases 
the gas ionization. Thus, the H+ and S"*" spatial dis- 
tributions are partially decoupled. The effects of ion- 
ization on the Ne{SB)/Ne{[Sll]) ratio are generally 
modest; they become important in the case of a radial 
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Fig. 3. True electron density structure in the four slices of NGC 1501 covered by the spectroscopic slit during "run 2" . 
The lowest density shown is A^e = 200 cm""^, while the highest one {Ne = 1410 cm~^) is reached by the approaching 
knot in PA= 30° (N-E sector), at an apparent distance of 11 arcsec from the central star. 



density profile which is multi-peaked, very asymmetric 

and strongly variable in direction; 
d: measurement uncertainty connected to the weakness of 
low ionization emissions. A symptom of the difficulties 
encountered by AUer and Epps (1976) comes from the 
non-detection of the [SII] red doublet at a second po- 
sition close to the nebular edge (25 arcsec from the 
nucleus, in PA= 223°), which is noticeably brighter 
than the first one. For this second position (interme- 
diate between our PA= 30° and PA= 55°) wc derive 
electron density peaks of 1050 — 1150 cm^'^ Note that 
these Ne{SB) values are only 10% higher than those 
obtained for the first position, indicating that the large 
flux observed in position 2 is mainly due to projection 
effects (i.e. thickness of the emitting layer). 

Most likely, to explain the Ne{SB)-Ne{[Sll]) discrep- 
ancy, we must look for a suitable combination of all the 
factors listed above. New, valuable information will come 
out of the detailed spatial reconstruction of the nebular 
structure, as performed in the next Section. 

5. 3-D spatial reconstruction 

5.1. Methodology 

Each tomographic map of NGC 1501 can be regarded as 
an image Ne{x, x'), where x is the distance from the plane 
of the sky crossing the central star and x' is measured 
along the slit. In order to recover the spatial electron den- 
sity distribution Ne{x,y, z), a suitable reference system 
is defined: we adopted x pointing toward the Observer, y 
toward East and z toward North. 

A data-cube of a given size N x N x N is initialized, 
filled with the raw data and re-mapped into cylindrical 
coordinates {x,r,a), where: 

r = ^/{y'^ + z"^); a arctan (^-^ (3) 

Once we have defined the pixel size p, we can fill up 
the n planes related to the spectrograms as: 

Neix,r,a)=NeJ-/-] . (4) 

Interpolation between the n observed planes is made in 

a linear fashion along the cylindrical coordinates. In other 

words, the data for a. > a, and a < ttj+i are given by: 

,r / N »T /a; x'\ da~ ,^ fx x'\ 6af 
Neix,r,a)=Ne.^-,-j ^+Ne,,. [-,-) ^(5) 

where: 

6a^=a - ai, Saj'=ai+i - a and Aa=ai+i - ai. 
Finally, one should consider the singularity my = z = 

0. These points are treated as an average of the neigh- 
bouring pixel-data for each plane defined by a given x. 



In the examples reported here we used N = 200, n = 8 
and p = 0.40 arcsec. 

5.2. Results 

A problem arises when presenting the A^e data-cube of 
NGC 1501: how can we render the 3-D on the page? The 
adopted solution is illustrated in the following figures. 

Fig. 4 shows the opaque reconstruction of the densest 
nebular regions (A''e > 900 cm~^) as seen from 12 direc- 
tions, separated by 15°. 

Fig. 5 is the same as Fig. 4, but for the weakest com- 
ponents {Ne > 300 cm~^); at even lower densities the 
roundish halo appears. 

Finally, Fig. 6 represents the optical appearance of the 
rebuilt-nebula seen from the same 12 directions of Figs. 4 
and 5. 

In these figures the line of view is identified by the 
zenith angle (6) and the azimuthal angle {ip), correspond- 
ing to a rotation through the first two Euler angles. Thus, 
for example, the upper-right image represents the rebuilt 
nebula as seen from the Earth (i. e. from (0,0)). 

Let's consider Figs. 4 and 5. Since the line of view of 
two adjacent images differs by 15°, each horizontal cou- 
ple constitutes a stereoscopic pair, allowing the reader to 
enjoy 12 surprising 3-D views of the nebula in as many 
directions (instructions are given in the caption of Fig. 4). 
The stereoscopic view is more difficult in Fig. 6, because 
of the diffuse, smooth nebular appearance. 

Hereafter we will call "direct stereoscopy" the 3-D view 
just described; in sub-section 6.3 an alternative method, 
"swapped stereoscopy" , will be introduced and compared 
with the "direct" one. 

For reasons of space, only twelve low resolution neb- 
ular structures (at two density cuts) arc presented, and 
as many projections. The complete data-cube can be di- 
rectly requested from the authors. Moreover, some movies 
of NGC 1501 (-1- other pleasantries) will be available soon 
from a dedicated WEB page. 

5.3. Reliability tests 

The first (and most important) consistency test of our 

reconstruction method is based on the comparison of the 
optical appearance of the true- and the rebuilt-nebula, 
when reduced to the same angular resolution. 

To be noticed that the spatial resolution of NGC 1501- 
rebuilt is not constant, but anti-correlated to the angular 
distance from the central star (this occurs because we ob- 
served the nebula at radially arranged position angles). 
Thus, taking into account the slit width, the number of 
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Fig. 4. Opaque reconstruction of the densest regions of NGC 1501 (iVe > 900 cm~'^) as seen from 12 directions 
separated by 15°. The line of view is identified by {0,^), where 6 is the zenith angle and -0 the azimuthal angle, 
representing a rotation through the first two Euler angles; the upper-right image is the rebuilt-nebula as seen from the 
Earth, i. e. from (0,0). Each horizontal couple constitutes a "direct" stereoscopic pair, allowing the reader to have 12 
3-D views of the nebula in as many directions (all together covering a straight angle). Instructions for beginners: to 
obtain the three-dimensional vision, look at a distant object and slowly insert the figure in the field of view (always 
maintaining your eyes parallel). Alternatively, you can use the two small dots in the upper part of the figure as follows: 
approach the page till the two dots merge (they appear out of focus); then recede very slowly (always maintaining the 
two dots superimposed) till the image appears in focus. A suggestion: be patient. The difficulties encountered in the 
beginning soon disappear and you will have exciting three-dimensional views of the nebula. 

Fig. 5. Same as Fig. 4, but for the weakest components (iVe > 300 cm^"^) of NGC 1501. To enjoy the stereoscopic 
view of the nebula in 12 directions, read the instructions contained in the caption of Fig. 4. 

Fig. 6. Optical appearance of NGC 1501-rebuilt from the same 12 directions of Figs. 4 and 5 (and at the same scale). 
Note that projection (O,-0)= projection(O,'0±18O°). The upper-right image corresponds to the rebuilt-nebula as seen 
from the Earth. In this case the stereoscopic view is quite difficult, because of the smooth, diffuse nebular appearance. 



position angles and the nebular angular extension, we have 
that the spatial resolution of NGC 1501-rebuilt varies from 
1.5 - 2.0 arcsec (at the central star position) to 3.0 - 3.5 
arcsec (at the nebular edge). 

The optical appearance of the true nebula under a see- 
ing of 2.0 and 3.5 arcsec is shown in Fig. 7. The comparison 
with NGC 1501-rebuilt (Fig. | upper-right image) is very 
satisfactory, both qualitatively and quantitatively. Indeed: 

1. all the macro- and micro-characteristics observed in 
NGC 1501 (general appearance, holes position, blobs 
and condensations distribution etc.) are faithfully re- 
produced by the model; 

2. the surface brightness profiles of the two nebulae agree 
to within ±10% along the position angles covered by 
the spectroscopic slit and to within ±20% at interme- 
diate directions. 

Other checks, eg. the comparison of the nebular fluxes 
at different radii and/or different positions, the iVe distri- 
bution etc. confirm that our technique adequately recovers 
the 3-D nebular shape. In conclusion. Figs. 4 to 6 closely 
represent both the structure and appearance of NGC 1501 
rotating around a N-S axis centred on the exciting star. 

6. Discussion 

In general, the foregoing spatial reconstruction performed 
in different emissions (i. e. for various ions) allows us to 
obtain the detailed 3-D ionization structure of a FN. This 
is particularly valuable for objects presenting complex 
morphologies (butterflies, poly-polars, multi-shells etc.) 
and/or large stratification effects and/or small scale con- 
densations (ansae, knots, FLIERS= fast, low ionization 
emitting regions, BRETS= bipolar, rotating, episodic jets 
etc.). 

Moreover, photo-ionization models at unprecedented 
resolution of the whole nebula and of the micro- and 



macro-structures can be produced by combining the true 
spatial structure of the expanding gas with the ultraviolet 
stellar flux. 

Finally, the 3-D reconstruction can strongly contribute 
in solving the PNe "Problem", i.e. the distance: only 
the detailed knowledge of both the dynamical and phys- 
ical nebular properties will give the right interpretation 
of the angular expansion measured in first and second 
epoch HST and/or radio imagery, thus obtaining precise 
expansion-parallax distances (Hajian et al., 1993, 1995; 
Kawamura & Masson, 1996; Terzian, 1997). 

From all these points of view NGC 1501 is a complete 
disappointment: it is quite faint (i. e. only the brightest 
emissions can be analysed) and homogeneous (no knots, 
jets or ansae) at high excitation and density bounded 
(stratification of the radiation is negligible). 

In spite of this, the 3-D analysis of our "ordinary" PN 
is important to introduce the wide potentialities of the 
spatial reconstruction in studying at small and large scales 
both the phenomenology and the physical processes con- 
nected with expanding nebulae (PNe, nova shells, super- 
nova remnants, nebulae around Population I Wolf-Rayet 
stars etc.). 

6.1. General 3-D structure of NGC 1501 

NGC 1501 is a high excitation, density bounded PN, 
and its spatial structure is the same in TVe, Af(H+) and 
Af(0++), other than a scaling factor. 

Low ionization regions (for example N+, identified by 
A6584 Aof [Nil]) are too weak for a detailed 3-D analysis. 
Probably they are localized in the external parts of the 
densest knots and condensations, where the shadowing of 
the ultraviolet stellar flux lowers the nebular ionization. 

The results obtained in PI (the main body of NGC 
1501 is a thin ellipsoid of moderate ellipticity deformed by 
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Fig. 7. Optical appearance of NGC 1501 (the true nebula) 
under a seeing of 2.0 arcsec (left) and 3.5 arcsec (right). We 
have first removed the stars from the TNG R frame, and 
later blurred the image by convolution with a Gaussian 
having FWHM=2.0 arcsec and 3.5 arcsec, respectively. 
The field and the orientation are as in Fig. 1. You can 
compare these images of the true nebula with the one of 
NGC 1501-rebuilt seen from (0,0), i.e. with the upper- 
right image of Fig. 6 (recall that the spatial resolution of 
NGC 1501-rebuilt is anti-correlated to the distance from 
the central star: 1.5 - 2.0 arcsec at the stellar position to 
3.0 - 3.5 arcsec at the nebular edge). 

Fig. 8. Optical appearance of NGC 1501-rebuilt seen from 
a) (25,75) i. e. along the major axis of the ellipsoid; b) from 
(65,-15) i. e. along the intermediate axis and c) from 
(—25, —15) i. e. along the minor axis. The thin ellipsoid's 
section and the axes contained in the plane of each figure 
are shown. 



a pair of large lobes along both the major and intermediate 
axes and by a multitude of smaller bumps spread on the 
whole nebular surface, making the general 3-D structure 
like a boiling, tetra-lobed shell) can now be improved and 
extended. 

The thin, inhomogeneous shell forming the main neb- 
ula (Fig. 4) is characterized by a dense, oblique belt defini- 
ing the "equatorial" plane, i. e. the plane containing the 
axes b and c of the central ellipsoid. The maximum den- 
sity (A^e ~ 1400 cm~^ for £( = 1) is reached in the direction 
of the minor axis. Only some arcs and knots belonging to 
the lobes appear in Fig. 4. 

A lumpy, quite homogeneous cocoon at low density 
completely embeds the main shell (Fig. 5). The com- 
plex radial structure of this envelope (sharp outwards and 
broad inwards) cannot be seen here and will be presented 
and discussed later (sub-section 6.4). 

The "boiling" structure of NGC 1501 is clearly visible 
in Figs. 4 and 5 for the shell's regions with higher spatial 
resolution, i. e. those projected on -or close to- the central 
star when the nebula is seen from (0,0). 

To highlight the "tetra-lobed" shape, in Fig. || we 
present the optical appearance of NGC 1501 as seen from: 

a) (25, 75), i. e. along the major axis of the ellipsoid; 

b) (65, —15), i. e. along the intermediate axis; 

c) (—25, —15), i. e. along the minor axis. 

In conclusion: the main body of NGC 1501 is a de- 
formed, almost oblate ellipsoid of moderate ellipticity 
(a ~ 44 arcsec, a/b ~ 1.02, a/c ~ 1.11), denser in the 
equatorial belt containing the axes b and c. Hydrodynam- 
ical processes (possibly Ray leigh- Taylor instabilities and 
wind interaction) partially swept up the lower density re- 
gions of the ellipsoid, generating both the hemispheric 



caps and the broad inward tail in the radial matter distri- 
bution. 

The ionized nebular mass {Mion, coinciding with the 
total nebular mass) is 0.15(±0.03) Mq and the dynamical 
age 5000(±500) years. 

As noticed in PI, a small, attached halo extending up 
to 34 arcsec from the central star surrounds the main neb- 
ula. The roundish, homogeneous appearance of the halo 
suggests that it represents photospheric material ejected 
by the FN progenitor in the late AGB evolution. More- 
over, the radial surface brightness profile, steeper than 
the r~^ law expected of a steady-state flow, indicates that 
the halo corresponds to a phase of "enhanced" stellar 
mass-loss. Finally, the electron density profile (Fig. 2 of 
this paper and Fig. 4 of PI) gives a mass-loss rate of 
3(±l)xl0^^ M0 yr^^ and an age of 1.0(±0.2)xl04 years 
(assuming Vexp = 20 km s^^; see Habing, 1996 and refer- 
ences therein). 

We claim that the halo represents the vestiges of the 
"superwind" phase (lasted 5(±l)xl0'^ years) which gen- 
erated the nebula. Most of the gas was later swept up 
by wind interaction, and was completely ionized by the 
UV flux of the evolving star, thus producing NGC 1501. 
Note that the mean mass-loss rate in the superwind phase 
resulting from the analysis of the main nebula, given by 
(Mjo„ /superwind duration), coincides with the value just 
obtained for the halo. 

At this point a quick comparison with NGC 40 (the 
first PN tomographically studied in detail, Sabbadin et 
al., 2000a) can be instructive. 

NGC 40 is an optically thick, very low excitation 
barrel-shaped nebula with thin arcs emerging at both ends 
of the major axis; it is powered by a luminous WC8 star 
presenting a large mass-loss rate. According to the "born- 
again scenario" proposed by Iben (1984; see also Blocker, 
1995 and references therein), the central stars of NGC 40 
and NGC 1501 suffered a final thermal pulse in the late 
post-AGB phase, ejecting the photospheric strata and ex- 
posing the C and O rich core. 

Moreover, both the evolutionary sequence [WC/ate] - 
[WCeaH^] - [WC-PG 1159] - [PG1159] ( |Hamann, 19"97 ; 
pK^oesterke et al., 1997 ) and the dynamical ages (3500 ±500 
yr for NGC 40 and 5000 ± 500 yr for NGC 1501) would 
indicate a earlier evolutionary phase for NGC 40. 

In 2000 ± 500 yr its size should be comparable with 
the present size of NGC 1501, but the structure will be 
quite different: a dense (optically thick?) equatorial torus 
+ extended, faint polar caps. 

Thus: NGC 40 will most likely evolve into a bipolar 
(butterfly?) PN, whereas NGC 1501 will keep its ellip- 
soidal structure, becoming fainter and fainter (in a few 
thousand years it will resemble A 43 and NGC 7094, 
two extended, filamentary, high ionization PNe excited by 
hybrid-PG1159 type central stars; see Rauch, 1999 and 
Feibelman, 2000). 
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The different evolution of tlie two nebulae can be ten- 
tatively explained by a more massive NGC 40 progenitor, 
which ejected a higher nebular mass with a larger density 
gradient between equatorial and polar regions, and created 
a higher mass central star (with respect to NGC 1501). 

Further support for this qualitative scenario comes 
from a preliminary analysis of a third "born-again" PN 
at an even earlier evolutionary phase, BD-f 303639 (the 
Campbell's star), observed with TNG + SARG in nine 
position angles at spatial and spectral resolutions of 0.7 
arcsec and A/AA ~ 115000, respectively. 

6.2. 2-D (morphology) - 3-D (structure) correlation 

Because of its fascinating variety and importance in most 
aspects of the advanced evolution of low and intermediate 
mass stars, the morphology of PNe has attracted many 
authors and several efforts have been made to search for 
a general classification scheme for the "morphological for- 
est" (Greig, 1972; Zuckerman & AUer, 1986; Bahck, 1987; 
Chu et al., 1987; Balick et al., 1992; StangheUini et al., 
1993; Corradi & Schwarz, 1995; Manchado et al., 1996). 

At the same time, a number of hydrodynamical and 
magnetohydrodynamical simulations were carried out and 
the theoretical nebulae compared with the observed ones 
(Mellema, 1997; Dwarkadas & Bahck, 1998; Garcia- Se- 
gura et al., 1999 and references therein). 

Our 3-D reconstruction method, allowing us to derive 
the detailed spatial structure of each object and to ob- 
serve it from all possible directions, revolutionizes the ap- 
proach to the morphological problem of PNe and opens 
new prospectives for understanding the 2-D (morphology) 
- 3-D (structure) relation. 

Let's consider Fig. 6, for instance, giving a represen- 
tative sample of the manifold morphologies of NGC 1501 
when changing the line of view. 

The optical appearance of the nebula mainly depends 
on the "latitude" . i. e. on the angle of view of the dense. 



In short: all the possible morphologies assumed by 
NGC 1501 when seen from different directions fall in the 
categories named C (centric PN) by Grcig (19"72|), Def (e l- 
liptical disk with filaments) by ^uckerman fc AUer (1986 ), 
E (elhptical) by [Bahck (1987D, ES (el liptical with inner 
filaments) by ^tangheUini ct al. (1994| ), E (elhptical) by 
Corradi & Schwarz (1995) and Es (elliptical with inner 
structure) by Manchado et al. (1996[ ) . They correspond to 
quite regular and symmetric shells, rings or disks prob- 
ably ejected by single, low mass progenitors (1.1 M© < 
Mms < 2.5 M©), excited by low mass central stars (0.56 
M© < Mcs < 0.60 M©) showing a galactic distribution 
and kinematics typical of the old Disk Population (Peim- 
bert & Serrano, 1980; Corradi & Schwarz, 1995; Maciel & 
Quireza, 1999 and references therein). 

In a certain sense, all this belongs to the past: the un- 
derstanding of the different mechanisms and physical pro- 
cesses forming and shaping PNe passes through the knowl- 
edge of the envelope structures. In their turn, these 3-D 
structures are usually presumed from a schematic clas- 
sification of the manifold nebular morphologies. Such a 
rough cascade process is overcome by our 3-D analysis, 
which directly provides the accurate spatial distribution 
(and ionization) of the expanding gas, thus removing any 
misleading camouflage due to projection. 

An example: while it is evident that in no case will 
NGC 1501 appear as a butterfly or a hourglass, the oppo- 
site could occur: when seen at particular directions (e. g. 
along -or close to- the axis of the lobes), a genuine butter- 
fly or hourglass PN (hke NGC 2346 or NGC 650-1) could 
mimic some NGC 1501 appearances (as suggested for K 
4-55 by Guerrero et al., 1996). 

This could also be the case for some famous and exten- 
sively studied PNe, like NGC 6720. The three-dimensional 
structure of the Ring Nebula has never been firmly de- 
fined, spanning from a complex toroid viewed close to the 
polar axis ( Minkowski fc Ostcrbrock, 1960| ; Louise, 1974 



Reay fc Worswick, 1977; Balick et al., 1992; Volk fc Leahy 



1993| ; [Bryce et al., 199| ; [Lame fc Poggc, 1994|) to an oblate 



inhomt pgeneous equatorial belt: at "low latitudes" (from 
(0.135 ) to f0.15) in Fig. 6) it is a broad disk prcscnt- 



ing inn er, amorphous structures. When seen almost pole- 



on (fr jm (0.60) to (0.105)). a sharp, oval ring appears. 1990|; fPascoli, 1990^ [Guerrero et al., 1997| ), to i flat ring 



bright er along the minor axis (corresponding to the pro- Pua fc Louise, 1970|) and to a hollow cylinder ( [Proisy 



or prolate spheroidal or ellipsoidal shell (Athcrton ct al 
19781 ; [PhUhps fc Reay, 1980[; [Kupferman 1983| ; [Masso^ 



iection lof the c axis of the ellipsoid). 



1974; Bachiller et al., 1989). No doubt the Ring Nebula 



A quick look through the main imagery catalogues 
([Acker et al., 1992[; [Schwarz et al., 1992| ; [Manchado et al 



1996| ; [Corny et al., 1999| ) indicates that NGC 1501-rebuih 



resembles M 3-30 when seen from (0,15), A 53 and NGC 
7094 from (0,30), IC 1454 from (0,45), A 73 from (0,60), 
A 70 from (0,105), NGC 4071 from (0,120) and NGC 6894 
from (0,135). 

We can add that, seen from other directions (not pre- 
sented here for reasons of space), NGC 1501 looks like 
several other nebulae contained in the above-mentioned 
catalogues. 



will be one of our next targets. 



6.3. Zooming in NGC 1501: local filling factor, ei, and 
density fluctuations in the main shell 

Figs. 4, 5 and 6 clarify that the big holes characterizing 
the optical appearance of NGC 1501 represent localized 
shell regions of lower density. This is even more evident in 
Fig. 9, showing the "swapped stereoscopic" nebular struc- 
ture for Ne > 700 cm"'^ as seen from (0,0), i.e. from the 
Earth. 
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Fig. 9. Opaque reconstruction of NGC 1501 for Ne>700 cm^'^, as seen from (0,-7; left) and (0,7; right), allowing the 
reader to have a "swapped stereoscopic vision" of the nebula from (0,0), i. e. from the Earth. To enjoy the 3-D view, 
cross your eyes till the two dots in the upper part of the figure perfectly coincide (in this way the left eye is forced 
to observe the right image, and viceversa). "Swapped stereoscopy" presents a great advantage with respect to the 
"direct" one: the image size can be as large as one wishes, thus achieving a better resolution than the "direct method" 
(where the centres of two adjacent images must be closer than the pupils separation). The drawback of "swapped 
stereoscopy" is the prohibitive physical size of a figure containing a large number of images (as in the case of Figs. 4 
and 5). Therefore, "direct stereoscopy" is convenient when presenting many pairs. Notice that in both methods a 
myope has an advantage, since he/she focuses closer. Caveat for beginners: when looking at a swapped stereoscopic 
pair for a long time, a slight trouble can arise (due to squint). Don't worry, but be careful. The field and the orientation 
are as in Fig. 1. 



The hole in PA= 70° at 11 arcsec from the central star 
corresponds to a "window" in the approaching nebular 
gas; through this "window" we observe the receding, dense 
material in the oblique belt representing the equatorial 
region of NGC 1501. The hole is clearly visible also in 
the following tomographic maps: Fig. 3 at PA= 80°, East 
sector and Fig. 5 of PI at PA= 55°, N-E sector. 

Since in the hole A'^e (receding gas) /Ae (approaching 
gas) > 1.6, we have flux(rec. gas)/flux(app. gas) > 2.5; 
i. e. to a first approximation we can neglect the contribu- 
tion of the approaching material and assume that we are 
observing only the receding one. In this way we have iso- 
lated a well-defined nebular portion, appropriate for study 
in detail at the maximum spatial resolution (TNG or, bet- 
ter, HST imagery). 

Figures 10a, b and c show the flux maps in the 15 x 15 
arcsec^ nebular region centred on the hole at 11 arcsec 
from the central star in PA= 70°, as given by our TNG 
broad-band R frame (Fig. 10a), by the same frame after a 
soft Lucy- Richardson restoration (Fig. 10b), and by HST 
(Fig. 10c). 

These figures stand out intensity fluctuations in the 
hole by a factor of 2.0, 2.5 and 3.0 for the TNG, TNG + 
Lucy-Richardson and HST images, respectively, indicating 
a clear correlation with the angular resolution. 

Let's take in the receding portion of the nebular shell 
observed through the hole: 

- the HST value as representative of the small scale flux 

variations, 

- the Ne radial profile derived from the 3-D reconstruc- 

tion, 

- the constancy of the electron temperature. 

Under these assumptions we obtain fluctuations of the 
density peaks of a factor 1.5(±0.2). 

Moreover, using a black-body spectrum for the exciting 
star and "normal" chemical abundances for the expanding 
gas (see Aller & Ckyzak, 1983), we have that: 

- the electron density peaks given by the [SII] red dou- 
blet are 1.10 - 1.20 larger than Ne{SB); 

- due to the incomplete recombination of S++ , the [SH] 
6717/6731 A density peaks are 10 to 30% lower than the 
true Ne values. 



Fig. 10. Flux distributions in the nebular region (15x15 
arcsec^) containing the "hole" at 11 arcsec from the star 
in P. A.= 70°, as obtained from our TNG R frame of NGC 
1501 (Fig. 10a), from our TNG R frame after the applica- 
tion of a soft Lucy-Richardson restoration (point spread 
function=:secing and ten iterations; Fig. 10b), and from 
the available HST imagery (Fig. 10c; we have assembled 
the four WFPC2 frames taken by H. Bond). The hole is 
the roundish region 8 arcsec in diameter at the centre of 
the 15x15 arcsec^ field. The field stars (truncated peaks) 
projected just outside the hole can be used as identifica- 
tion markers in Fig. 1. 

This is valid for a wide range of nebular distances 
(1.0 < D[Kpc] < 2.0) and stellar fiuxes (3.0 < 
logL»(L0) < 3.6; 4.90 < logT, < 5.10). 

Very similar results come from the analysis of other 
"holes" present in NGC 1501, for instance the one at 15 
arcsec from the star in PA= 125°, corresponding to a 
"window" in the approaching gas, and the one at 11 arc- 
sec from the star in PA= —13°, representing a "window" 
in the receding material (this last hole is visible in Fig. 1 
and in Fig. 3 at PA= 170°, Northern sector). 

In spite of the heavy assumptions, the foregoing re- 
sults would indicate small-scale density fluctuations (com- 
bined with ionization) as the main cause of the A^e([SII]) — 
Ne{SB) discrepancy reported in Section 4. 

This implies that: 

- the local filling factor in the main shell is ei = 

0.50(±0.1), i. e. lower than the value derived from 
{Ne{SB)/Ne[SIl]f = 0.66(±0.1); 

- our surface brightness electron densities should be in- 

creased by a factor 1.4(±0.2) to obtain the density 
peaks of the emitting gas in the main shell. 

It is evident that new insights into the variation on a 
small scale of the physical conditions are needed to deepen 
our understanding of the above-mentioned questions; this 
mainly concerns the electron temperature, whose changes 
affect in different ways the emissivities of recombination 
and forbidden lines (a lower Te favouring recombinations; 
see Peimbert, 1994, for a recent review). 
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As suggested in PI, an even lower value of e; is ex- 
pected in the inwards tail of the radial Ne distribution, 



by the large value of M (logAf = -6.28 MQ yr~^) ob- 
tained for the WC4 nucleus of NGC 1501 by |Koestcrke~et 



this knatty and inhomogcncous region arising from hydro- al. (1997) using the standard atmosphere model of Wolf- 



dynamical processes. 

Moreover, the application of eq. (1) to the whole neb- 
ula gives a global filling factor e < 0.25 (corresponding to 
the parameter normally used to estimate the fraction of 
the total nebular volume filled by the ionized gas; see Bofh 
& Stanghellini, 1994, and references therein). 

The zooming capability implied in our 3-D analysis is 
a powerful tool to overcome the existing gap between the 
potential resolution of the modern photo-ionization mod- 
els and the actual resolution when applied to real nebulae. 
The first (for example CLOUDY; Ferland et al., 1998) can 
be as detailed as one wishes simply using proper input 
parameters (stellar flux distribution, nebular dimensions, 
density law etc.), whereas their practical application is 
limited by observational and reduction procedures and by 
projection effects: up to now, all efforts intended to iso- 
late a well-defined nebular portion (long-slit spectra anal- 
ysed at different positions, short-slit echellograms etc.; see 
Barker, 1991, Perinotto & Corradi, 1998 and Hyung et al., 
2000), in the best cases give emission line intensities which 
are a mixing of both the approaching gas and receding gas 
spectra. 

The zooming yield is correlated to the spatial resolu- 
tion of the rebuilt nebula (i. e. to the spatial and spec- 
tral resolutions of the echellograms) ; it is quite modest in 
the present study of NGC 1501, performed at spatial and 
spectral resolutions of 1.5 — 2.5 arcsec and 22000 — 25000, 
respectively. The output will increase when analysing the 
symmetric knots of the quadrupolar PN IC 4634, the ansae 
and the FLIERS of NGC 7009 (the "Saturn") and the 
"wings" of the butterfly PN HE 5 (observed with ESQ 
NTT -I- EMMI at spatial and spectral resolutions of 1.0 
arcsec and 60000, respectively). 

6.4^. Radial electron density profile and open questions 

The intersection of the rebuilt-nebula, seen from (0,0), 
with the plane perpendicular to the line of sight and cross- 
ing the central star, gives the Ne map shown in Figure 11, 
representing the radial electron density profile (for e; = 1) 
extended to all position angles. 

Fig. 11 shows in detail the most intriguing character- 
istic of the gas distribution in NGC 1501, i. e. the broad 
inwards tail present in the radial density profile, probably 
arising from hydrodynamical processes (Rayleigh- Taylor 
instabilities and winds interaction; Capriotti, 1973; Kahn 
& Breitschwerdt, 1990, Garcia-Segura et al., 1999). 

The first process worked for about 1000 yr in the 
early evolutionary phases, when the nebula was ionization 
bounded (up to R ~ 0.06 pc). 

The second one, being independent of the optical thick- 
ness, acted for much longer (for most of the nebula's life, 
i. e. 5000 yr) and is probably still present, as suggested 



Rayet stars. 

This simplified scenario is complicated by the peculiar- 
ities of the powering nucleus of NGC 1501 (very hot W-R, 
rich in C and O, hydrogen deficient, pulsating), suggest- 
ing that the star suffered a thermal pulse in the late post- 
AGB evolution. Its path in the HR diagram is quite uncer- 
tain, given the large number of parameters involved: phase 
in the nuclear burning cycle at which the star departs 
from the AGB; mass of the residual hydrogen-rich stellar 
envelope; mass-loss rates during the quiescent hydrogen- 
burning phase, the high-luminosity phase following the he- 
lium shell flash and the low-luminosity quiescent helium 
burning phase; efficiency of convective overshoot in mixing 
hydrogen-free material out to the stellar surface; chemical 
mixing due to stellar rotation (Schonberner, 1979; Iben, 
1984; Blocker, 1995 and references therein). 

Thus, the detailed analysis of the mechanisms and 
physical processes forming and shaping the nebula is de- 
ferred to a deeper knowledge of both the star evolutionary 
parameters and the gas structure. 

A final remark concerns the rarity of the radial den- 
sity profile shown by NGC 1501: no evidence of an ex- 
tended inward tail is present in the other two dozen PNe 
up to now spectroscopically observed by us with the Asi- 
ago 1.82m telescope + Echelle, ESQ NTT + EMMI and 
TNG + SARG. 

A good candidate could be NGC 6751; the complex 
morphology of this nebula (excited by a WC4 star) was 
shown by Corradi et al. (199^) and spectrosco pically stud- 
ied at high dispersion by Chu et al. (1991), who iden- 
tified six different components having as many spatio- 
kinematical characteristics, but didn't give information on 
the radial density profile of the ionized gas. 

All this, on the one hand confirms and extends the 
results already obtained in PI, and highlights the poten- 
tialities of our reconstruction method, on the other hand 
leaves unresolved many exciting questions, such as: 

- what is the small-scale matter distribution in the 
main shell? And that in the inward tail? 

- did the hydrogen-depleted stellar wind create chem- 
ical composition gradients across the nebula? 

- how large are the electron temperature and turbu- 
lence fiuctuations in the shell? 

- what are the physical conditions, dynamics, ionic and 
total abundances in the outer, roundish halo? 

- is there any observational evidence supporting the 
"born-again scenario" proposed by Iben et al. (1983; see 
also Blocker, 1995, and Herwig et al., 1999) for nonradial 
g-mode pulsators? 

Most of these (and many other) questions can be an- 
swered by a deeper spectroscopic study at higher spatial 
and spectral resolutions (and at more position angles). 
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Fig. 11. Ne map and isophotal contours in the sky plane 
crossing the central star when the nebula is seen from 
(0,0), showing the pecuhar radial density profile of NGC 
1501. The density contours cover the range 200 cm""^ 
to 1200 cm"^ and are spaced by 100 cm~'^. This figure 
corresponds to the zero-velocity pixel column extended 
to all position angles (when using high resolution, slit 
spectroscopy), and to the rest-frame (for Fabry- Perot im- 
agery). 



6.5. Application of the 3-D analysis to expanding nebulae 
in general 

PNe, nova shells, nebulae ejected by symbiotic stars, bub- 
bles surrounding early spectral type main sequence stars, 
supernova remnants (SNRs), shells around Population I 
W-R stars all consist of expanding masses of ionized gas, 
and the many facets of the 3-D analysis can be success- 
fully applied to them for studying morphology, dynamics, 
physical conditions, photoionization model and evolution- 
ary status. 

In all cases the starting point is high resolution slit 
spectroscopy (or Fabry-Perot interfcrometry) . Given the 
large range of expansion velocities shown by the different 
classes of nebulae (a few dozen km s^^ for PNe to several 
hundred km s~^ for young SNRs and nova shells), the 
term "high resolution" is here intended relative to 'Vexp- 

Thus, the first parameter to be considered is the "rel- 
ative resolution", RR=Vexp/AV, AV being the spectral 
resolution. 

The best RR value is always a compromise between 
two contrasting factors, dispersion and detection: a bright 
nebula can be observed at large RR in both strong and 
weak lines, whereas only the brightest emissions are de- 
tectable in a low surface brightness object (even at small 
RR). 

An indicative lower limit to the applicability of the 3- 
D analysis can be put at RR=3. This means that slow 
expanding nebulae (like PNe) need high dispersion spec- 
tra (A/AA> 20000), whereas low dispersion spectroscopy 
(A/AA ~ 1000) is preferable for the fastest objects. The 
latter is clearly illustrated by Lawrence et al. (1995), who 
derived the three dimensional model of the emitting gas 
in the Crab Nebula, in Cassiopea A and in the remnant 
of GK Per (Nova Persei 1901) by means of low resolution 
(200 < A/AA < 1500) Fabry-Perot imagery. 

A second parameter to be considered in the 3-D anal- 
ysis is the spatial resolution along the slit (Ad). Also in 
this case: 

- it is intended "relative to the angular dimensions d 
of the nebula" ; 

- the larger the "relative spatial resolution" SS = 
d/Ad, the better the reconstruction; 

- as a rule of thumb, SS must be larger than 3. 



A detailed discussion of the various contributions given 
by the 3-D analysis in solving the manifold problems con- 
nected with the different classes of expanding nebulae is 
beyond the aims of this paper and is left to the educa- 
tion, experience and imagination of the reader. In fact, 
our intent is exactly this: to excite the interest of some 
other people, exposing them to tackle the many problems 
and the wide applications of the 3-D methodology (we 
feel quite alone at present... "rari nantes in gurgite vasto" 
paraphrasing Virgil, Aeneid I, 118). 

7. Conclusions 

This paper introduces a general reduction procedure giv- 
ing the detailed spatial structure of the ionized gas in ex- 
panding nebulae. 

Its application to NGC 1501 indicates that this high 
excitation PN is an almost oblate ellipsoid of moderate 
ellipticity (a ~ 44 arcsec, a/b ~ 1.02, a/c ~ 1.11), denser 
in the equatorial belt containing the axes b and c. A lumpy, 
quite homogeneous cocoon completely embeds the main 
shell. 

A series of images is presented, showing stereoscopic 
views of the 3-D gas structure and the nebular appearance 
when changing the line of sight. 

The morphology-structure correlation is discussed, the 
small scale density variations in the main shell and the 
peculiarities of the matter radial profile are analysed. 

NGC 1501 is a quite faint, high ionization, density 
bounded PN, its 3-D study is simple and straightforward 
(at our spatial and spectral resolutions). 

We are now involved in a much more ambitious project: 
the spatial mapping in different ions (H+, He+, He"'"+, 
0°, 0+, 0++, N+, S+, Ne++, Ar++, Ar+++ etc.) and the 
photo-ionization modelling of a dozen bright PNe in both 
hemispheres; each target observed at several (at least 9) 
position angles with ESQ NTT + EMMI (spatial resolu- 
tion = 1.0 arcsec, spectral resolution = 60000) and TNG + 
SARG (spatial resolution = 0.7 arcsec, spectral resolution 
= 115000). 
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